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Abstract 

Most models of Bose-Einstein correlations in multiple particle produc- 
tion processes can be ascribed to one of the following three broad classes: 
models based on the original idea of the Goldhabers, Lee and Pais, hydro- 
dynamic models and string models. We present for discussion some basic 
questions concerning each of these classes of models. 

1 Introduction 

Since the pioneering work of G. Goldhaber, S. Goldhaber, W. Lee and A. Pais, 
published over forty years ago Q and known as the GGLP model, Bose-Einstein 
correlations in multiple particle production processes have been studied in hun- 
dreds of papers. Many references can be found in the recent review articles 
[| and g. These correlations have been popular in particular for two reasons. 
They give impressive bumps in two-particle and many particle distributions 
and, if the regime of Einstein's condensation can be reached, there are even 
more spectacular phenomena waiting to be discovered [Q , jsj , 1^ , Q . What is 
more, Bose-Einstein correlations are believed to yield important information, 
which seems hard, if not impossible, to obtain by other means. They have been 
used to find the sizes and shapes of the regions where hadrons are produced, as 
well as to obtain detailed information about the evolution in time of the hadron 
production processes (see the reviews [D, |^). 

There is little, doubt, however, that the problem is hard. As an example of 
a sceptical opinion let me quote one of the creators of this field of research, G. 
Goldhaber, who said at the Marburg conference in 1990: "What is clear is that 
we have been working on this effect for thirty years. What is not as clear is 
that we have come much closer to a precise understanding of the effect" . Every- 
body agrees that the GGLP paper was very important and various extensions 
of the model proposed there are still being used. There are, however, other 
approaches. The GGLP model contains static sources of particles. In the more 
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recent "hydrodynamic" models the flow of the sources is of great importance. 
Another very promising approach are string models, where the random phase 
assumption used in previous models is not necessary and the description looks 
closer to QCD. 

In the present paper we will characterize the GGLP models, the hydrody- 
namic models and the string models, stressing in each case open and potentially 
important problems. 



2 GGLP models 

Let us consider two identical bosons, e.g. two 7r+ mesons, created: one at point 
fi and the other at point r2- If the two bosons were distinguishable, a crude 
approximation for the probability amplitude of observing both of them at point 
r could be 

j^^ _ gi<^i+ipi-(r-fi)gj<^i+ip2-(f-f2) ^-[^-j 

Here the interaction between the two bosons is neglected, so that the two-particle 

amplitude is a product of single particle amplitudes. Only the phase factors are 
kept and each phase is the sum of the phase obtained by the boson at birth and 
of the phase acquired while propagating with given momentum from the birth 
point to point r. 

For identical bosons, however, this amplitude has not the right symmetry 
with respect to the exchange of the two bosons and the least one must do is to 
symmetrize it. Thus for identical bosons the corresponding approximation is 

_ ^_gj(0l+<?i2)+i(pi+P2)-r ^g-j(pi-fi+p2-f2) _|_ g-i(P2-fi+pi-f2)~j ^2) 

The probability distribution for momenta is proportional to 

|Ap = l + cos((pi-p2)-(ri-f2)). (3) 

In order to make use of this expression it is necessary to average it over the non 
measured production points fi,f2- In the GGLP paper the averaging was over 
the space distribution of sources p{ri;R)p{r2;R), where i? is a parameter with 
dimension length, which was interpreted as the radius of the production region. 
Various generalizations, modifications and extensions followed, but let us use 
this simple variant to make some general remarks. 

In the GGLP model the distribution for pairs of particle momenta depends 
only on the momentum difference q = pi — P2- This is in violent contradiction 
with the data, but GGLP found a clever way out. The distribution for the 
unsymmetrized amplitude is flat. Therefore, the result can be just as well 
interpreted as a prediction for the ratio of the actual momentum distribution 
to the distribution for distinguishable bosons. Further we denote this ratio by 



R{PX,P2) = -> x|2 (4) 



I^07i.j72j|- 
|Ar,(pi,p2)|^ 

Now, the momentum distribution is not assumed to be independent of the sum 
of the momenta. It is enough to make the much weaker and more reasonable 
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assumption that the dependence on this sum can be factored out and cancels in 
the ratio. A weh known difficulty with this approach is that the distribution for 
distinguishable tt^-s, say, cannot be obtained from experimental data without 
further assumptions. GGLP assumed that the distribution for tt+tt" pairs can 
be used instead. There have been many other proposals (cf. e.g. Q and 
references contained there), but none is fully satisfactory. 

For any nonsingular averaging process the average cosine must be close to 
one for |(7| « and very small for large values of Therefore, the ratio 
R(pi,P2) decreases, though not necessarily monotonically, from values close to 
two for small values of \q\ to values close to one for large values of jifl. This 
gives the characteristic bump in R{pi,p2) for small values of g^. If i? is the 
only dimensional parameter available, the width of this bump must, for simple 
dimensional reasons, be proportional to R^^. Thus, the main qualitative re- 
sults of GGLP are much more general than their specific choices of the weight 
functions p{r;R). Nevertheless, they are not quite general. 

It is well known from optics that, whether photons bunch or not, depends 
on the type of source and not only on the fact that they are bosons. Photons 
can antibunch just as well. In order to illustrate this point within the GGLP 
type models, let us assume that the amplitude A has an additional factor, which 
equals one, if the product (ri — r2)(pi — P2) > and minus one otherwise. This 
factor changes sign, when the momenta of the two bosons are exchanged. There- 
fore, the squared modulus of the properly symmetrized production amplitude 
is 

|A|2 = l-cos((pi-p2)-(ri-r2)). (5) 

and we get a hole instead of the bump in the small g*^ region. Admittedly this 
model is not realistic. Its purpose is only to indicate a possibility. This may 
be interesting in view of the LEP results concerning Bose-Einstein correlations 
in e~^e~ annihilations, where two W bosons are simultaneously produced. It 
seems that identical pions originating from the decay of a single W exhibit the 
usual bump attributed to Bose-Einstein correlations, while these correlations 
are absent, or very weak, for pairs of identical pions, when each pion originates 
from a different 1^ §, @, In the GGLP model the bump results from 

the assumptions that pion pairs produced in different pairs of points add inco- 
herently. Mild modifications of this assumption Q , can affect the size of the 
bump, but do not eliminate it. It would be interesting to check, whether the 
GGLP assumptions could be modified so as to predict the bump for some, but 
not for all pairs of identical mesons produced in a multiple particle production 
event. 

3 Hydrodynamic models 

It is not possible to express the ratio R{p,p') in terms of the single particle 
momentum distributions. In the hydrodynamic models, as well as in GGLP 
models, one makes, however, an assumption, which makes it possible to ex- 
press this ratio in terms of the diagonal and off-diagonal terms of the single 
particle density matrix in the momentum representation. It is convenient to 
formulate the hydrodynamic models in terms of source fimctions S{X,K). The 
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source function ||l2|, is related to the single particle density matrix in the 
momentum representation by the formula 



In this formula 



K 



1 



(p + p); 



P- 



■P ■ 



(6) 



(7) 



Here K,q,p,p' are fourvectors, but in order to calculate the density matrix, we 
need only their values corresponding to the momenta p, p' being on their mass 
shells. X is an integration variable, which is associated with the position of the 
sources in space-time. The physical interpretation of X may be helpful, when 
trying to find the source function. It is irrelevant for the calculation of the 
density matrix, once the source function is known. 

There is an infinity of different source functions, which all give the same 
density matrix and consequently the same predictions for the ratio R{p,p'). For 
instance, one could put 

S{X,K)^W{X,K)d{Xo), (8) 
where W{X, K) is the well- know Wigner function satisfying the relation 



P{P,P) 



'iq-X 



W{X,K)(fX 



and 



X=-{x + x'). 



(9) 



(10) 



This source function gives the correct density matrix by construction, but it 
corresponds to a most unlikely scenario, where all the particles are created si- 
multaneously at Xq — 0. Since our aim is to find the correct density matrix, this 
source function would be fine, in spite of the unlikely physical picture attached 
to it. The problem is, however, that finding the Wigner function is not any 
easier than finding the density matrix in the momentum representation. The 
hope is that using a source function, which corresponds to a plausible scenario 
for the production process, we will be able to use more efficiently what we know 
about particle production in order to find the source function. 

As an example of this approach let us consider the model reviewed in 
The source function is postulated in the form 



S{X, K) = CrriT cosh(?/ — i]) exp 



rriT cosh y cosh r/f — Tj^^XKt sinh T]t 



exp 



T 



2i?2 2(A77)2 2(Ar)2 



(11) 



In this formula 
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Vt = Vf^; T = - z2; 

m"^ ^ m'^ + K'^ ; = (12) 

z is parallel to x\\ and x is parallel to K. The source function depends on 
six free parameters (R,T,rif, Atj^tq, At) and, moreover, contains the normal- 
ization constant C, but each piece of the source function has a clear physical 
interpretation. Therefore, fixing these parameters from the data yields directly 
interesting physical information. We will illustrate this important point using 
a fit to the NA49 data on Pb-Pb scattering at 158 GeV/c per nucleon The 
parameter R is the transverse radius of the tube, from which the final hadrons 
are emitted. The result R « 7fm is about twice the radius resulting from from 
the known radii of the lead nuclei. This is evidence for a significant trans- 
verse expansion, before most of the hadrons are produced. The parameter rjf 
governs the transverse rapidity of the sources. The value obtained rjf 0.35 
corresponds to transverse velocities reaching the velocity of sound in the plasma 
(1/3), which looks very reasonable. The parameter T occurs as the temperature 
in a Boltzmann type factor. Its fitted value T f« 130fm is significantly lower 
than the temperatures obtained, when fitting the chemical composition of the 
final state hadrons, which indicates that during the expansion the stuff cools 
down. Another interesting comparison is that of the fitted values tq « 9fm 
and At « 1.5fm. The parameter tq is the typical time between the moment of 
collision and the moment, when a hadron is produced. The parameter At is the 
duration of the time, when hadrons are produced. The fact that tq ^ At means 
that all the hadrons are produced in a short time interval after a relatively long 
incubation time. Unfortunately, as stated by the authors the parameter At 
is poorly constrained by the data, so that this conclusion is not as solid as the 
others. 

As seen from this example, given a model one can obtain from the data 
much important information. An open problem is, however, how stable are 
these conclusions when models change. As seen from the expression of the 
source function in terms of the Wigner function, one can fit perfectly the data 
assuming that all the particles are produced exactly simultaneously. It is just 
as easy to get a perfect fit assuming that the particles are produced only on 
the surface of a sphere, or only on the surface of a cube. These alternative 
models are so implausible physically that there is little doubt they should be 
discarded. The question is, however, how many physically plausible models can 
fit the data, while giving completely different descriptions of the hadronization 
process? 

4 String picture 

The string model for Bosc-Einstein correlations , , [|l6j has not yet been 
developed to the point, where it could be compared quantitatively with the data. 
It is, however, much more ambitious than the models described above. Instead 
of phenomenological assumptions about sources and incoherence, it gives a well 
defined amplitude for the production of particles with momenta pi, . . . This 
amplitude is a plausible approximation to QCD. We will consider only the 1-1-1 
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dimensional version of the string model, which seems to contain all the main 
ingredients of this approach, while it is much simpler than the 3+1 dimensional 
version. In fact the 1+1 dimensional model has been recently analytically diag- 
onalized ||l^, though in the version without the Bose- Einstein correlations. 

Let us consider a final state consisting of hadrons with momenta pi, . . .pn- 
To this final state the model ascribes a polygon in the (z, t) plane. The sides 
of this polygon are the trajectories of the various partons existing between the 
moment of e~^e~ annihilation and the moments, when the hadrons are formed. 
The partons are considered massless and moving with the velocity of light, 
therefore the sides of the the polygon form angles ±45° with the t and z axes. 
Let us put the e~^e~ annihilation point at the origin of the coordinate system. 
At this point two partons, a quark and an antiquark, are formed flying along 
the z axis, away from each other. Their trajectories form the first two sides of 
the polygon, both starting at the origin, one going to the right and upwards, 
the other going to the left and upwards. The partons are end points of a colour- 
string. Thus the string sweeps the surface of the polygon. The energy of the 
string E is connected to its length L by the formula E = kL, where k is a 
constant known as the string tension. Thus, while the quark and the antiquark 
fly away from each other and the string expands, there is a force reducing the 
energy of the two partons and finally the directions of their motions get reversed 
starting another pair of the sides of the polygon. In the meantime the string at 
any point between the endpoints can break producing a quark and an antiquark, 
which form another pair of sides of the polygon. Since any segment of the string 
has a quark at one end and an antiquark at the other, it is easily checked that, 
except for the original two partons, all the quarks fly in one direction and all 
the antiquarks in the other. From time to time a quark meets and antiquark. 
Then the two form a hadron with two-momentum equal to the sum of the two- 
momenta of the two meeting partons. Thus the polygon contains the following 
elements: the vertex at the origin, two turning points of the original partons, n 
vertices, where the n hadrons were formed and n — 1 vertices, where the string 
broke. One finds that the two-momentum of a hadron is determined by and 
determines the lengths of the two sides of the polygon adjacent to the vertex 
where the hadron was formed. 

The probability amplitude for producing the state {pi, . . . ,pn) depends on 
the area A of the polygon 

M(pi,...,p„)^e'«^. (13) 

The imaginary part of ^, ib, gives the probability distribution well-known from 
the LUND model 

\M{p,,...,p,,)\' ^e-"^. (14) 

The imaginary part, believed to be close to the string tension k, is important 
for the description of the Bose-Einstein correlations. In order to describe these 
Bose-Einstein correlations the amplitude ( p^ is symmetrized very much like in 
the GGLP approach and a qualitatively satisfactory description of the correla- 
tions is obtained. Since, however, the amplitude being symmetrized is not the 
GGLP one, there are some significant new points. 

Symmetrization means summing over all the permutations of identical hadrons. 
Let us concentrate on an exchange of two 7r+ mesons. Each of them is produced 



6 



at some hadronization vertex of the polygon. In order to perform the exchange, 
one has to cut off the two pionic vertices together with their adjacent sides of 
the polygon and to glue them back, pion one in the position of pion two and 
pion two in the position of pion one, so as to obtain again a closed polygon. The 
new polygon has in general a different area A' ^ A. If the new area is much 
larger than the area before the exchange, the contribution from the interference 
with the permuted amplitude is negligible. One reason, familiar from the GGLP 
model, is that the relative phase is a rapidly varying function of momentum. 
The new fact is, however, that the modulus of the permuted amplitude is addi- 
tionally suppressed by the factor e~^'^^~^ In order to obtain small changes 
of the area, it is advantageous to exchange pious, which are close to each other 
counting along the perimeter of the polygon, or equivalently, which have similar 
momenta. This is the reason for the familiar bump for pi p2 

Let us quote two interesting qualitative predictions of this model. There 
should be a difference between the Bose-Einstein correlations for pairs of charged 
pions, say tt+tt"'', and corresponding correlations for pairs of neutral pions. The 
reason is that two tt^-s can be formed at two adjacent hadronisation vertices of 
the polygon, while two 7r+-s cannot. There should also be a difference between 
the Bose-Einstein correlations for pairs of mesons originating from the same 
string and pairs of mesons originating from different strings. For the latter 
situation the present model is clearly not applicable. It has been suggested (Q 
and references quoted there) that perhaps there are no correlations between 
mesons from different strings, which would explain the observations of Bose- 
Einstein correlations for pions from decays of pairs of W bosons. 

An interesting question raised by Bowler is the relation of the string model 
to the GGLP model. Bowler proposed Q a model closely related to the GGLP 
model, which looks very similar to the string model and, according to Bowler, 
gives also very similar predictions. The model used by Bowler contains a dis- 
tribution of sources, which depends not only on space-time points, as in the 
GGLP paper, but also on the momenta of the produced particles. Such models 
have become popular after the work Yano and Koonin and are sometimes 
called Yano Koonin models. It is not clear, however, what constraints must 
be imposed on the distributions, in order to make them consistent with quan- 
tum mechanics. Therefore, occasionally this approach gives inconsistent results 
l ^lf . Bowler's model deviates in two ways from the string model. It yields the 
probability amplitude as an exponential in the area, but this is a modified area, 
where certain regions are counted more than once. Moreover, in both models, 
in order to get the inclusive fc-particle distribution it is necessary to integrate 
out the momenta of the remaining particles, but in Bowler's model the integra- 
tion region is different from that used in the string model. According to some 
unpublished numerical calculations by Bowler the two deviations nearly cancel. 
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